o 
o 

(N 
X) 

tin 

in 

(N 



o 
o 

u • 

Oh- 

a 

o 
o 



> 
o 

00 

in 

(N 

o 
o 

S 

O 

o 



Adiabatic quantum computation with flux qubits, 

first experimental results. 

S.H.W. van der Ploeg, A. Izmalkov, M. Grajcar, U. Hiibner, S. Linzen, S. Uchaikin, Th. Wagner, A. Yu. Smimov, 
A. Maassen van den Brink, M.H.S. Amin,A.M. Zagoskin, E. Il'ichev and H.-G. Meyer 

(Invited Paper) 



Abstract — Controllable adiabatic evolution of a multi-qubit 
system can be used for adiabatic quantum computation (AQC). 
This evolution ends at a configuration where the Hamiltonian of 
the system encodes the solution of the problem to be solved. As 
a first steps towards realization of AQC we have investigated 
two, three and four flux qubit systems. These systems were 
characterized by making use of a radio-frequency method. We 
designed two-qubit systems with coupling energies up to several 
kelvins. For the three-flux- qubit systems we determined the 
complete ground-state flux diagram in the three dimensional flux 
space around the qubits common degeneracy point. We show that 
the system's Hamiltonian can be completely reconstructed from 
our measurements. Our concept for the implementation of AQC, 
by making use of flux qubits, is discussed. 



I. Introduction 

THE concept of quantum computation by adiabatic evo- 
lution, i.e. Adiabatic Quantum Computation (AQC), was 
introduced by Farhi and coauthors [1]. It is an alternative to the 
"standard" concept of quantum computation where one has to 
construct a universal set of gates with long coherence times. As 
the latter requirement is not easy to achieve in superconducting 
quantum bits (qubits), presently the longest times achieved are 
of the order of several /is and only achieved at an optimal 
point [2], [3], AQC seems to be better as the system is only 
required to stay in the ground state. This led Kaminsky and 
co-authors to propose a scalable superconducting architecture 
for adiabatic quantum computation [4], where they make use 
of 3 Josephson Junction persistent current qubits [5]. 

In the AQC concept the solution of some non-polynomially 
(NP) hard problem is encoded into the ground state of a 
complex multi-qubit Hamiltonian Hp . The parameters of such 
a Hamiltonian can be controlled continuously so in this sense 
this method is similar to analog computation. The computation 
itself contains three majors steps: (A) preparation of the system 
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Fig. 1. Experimental setup for the demonstration of adiabatic quantum 
computation in a three-qubit-system, see section |v] The qubits are coupled to 
a resonant tank circuit for the readout. The flux bias to the qubits is provided 
by a dc current IhT applied to the coil Lt and bias currents Ibi,2 applied to 
additional bias wires, allowing full control over the individual qubits fluxes 
fi. In the case of more qubits, more current sources would be used. 



in a well-known state with Hamiltonian H(t = 0) = Hi at 
time t = 0, (B) adiabatic evolution during the time interval 
^caic towards the problem Hamiltonian Hp and (C) readout 
of the ground state of the problem Hamiltonian, which gives 
the answer. Therefore, in order to build a quantum computer 
based on AQC one has to be able to construct a qubit system 
with a controllable Hamiltonian and one has to be able to read 
out this system once the problem Hamiltonian Hp is encoded 
into it. 

There has been an experimental implementation of an 
AQC algorithm in an NMR three-qubit quantum computer by 
Steffen et al [6]. By using their computer they solve the so- 
called MAXCUT problem. This problem is one of the Non- 
Poly nomially hard problems of Graph theory where the object 
of the problem is to find the maximal "cut" of a graph [7] 
and is mathematically equivalent to finding the ground state 
of some Ising-type Hamiltonian. Some of the authors have 
proposed how to implement and readout such an algorithm 
making use of flux qubits in a setup like the one shown in 

Figmm. 

In this paper we review our recently obtained results from 
coupled flux qubits using our and implemented method for 
their characterization. Also we will show some of our results 
on multi qubit systems. 

n. Coupling of qubits 

In order to implement any protocol for quantum informa- 
tion processing, qubits should be coupled. Three Josephson 
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junction persistent current flux qubits [5] can be coupled 
inductively just by placing them next to each other, however 
the achievable coupling strength is rather weak [9]. A larger 
coupling strength can be achieved by making use of the 
kinetic inductance of a common leg between them [10]. Even 
stronger coupling can be achieved by inserting a Josephson 
junction in this common leg [11]. The latter also provides 
the possibility of ferromagnetic coupling, whereas the other 
methods only provide anti-ferromagnetic coupling [12]. By 
insertion of an additional coupler loop between the qubits, 
one could realize controllable/switchable coupling between 
them [13]. This is particulary important for the realization of 
scalable superconducting architecture for AQC. Controllable 
coupling enables one to construct a Hamiltonian where both 
the energy bias and the coupling terms Jij can be set in-situ. 

A coupled A^-flux-qubit system can be described by the 
Hamiltonian: 

N 

H, = - Y}^iaf + AicrW] + ^ Jaofai^^ (1) 

i i<j 

Here is the tunneling amplitude, or half energy splitting, 
of qubit i at /i = and Si = ^oIp,ifi gives the energy 
bias applied to qubit i in terms of its persistent current Ip^i 
and the normalized flux bias fi = ^x,i/^o — 1/2 (<l>o is the 
flux quantum). The symbols crx\(Jy\ai^'^ denote the Pauli 
matrices of the i-th qubit. The coupling energy is given by: 

where Mij represents both the magnetic inductance between 
qubits i and j, as well as the kinetic inductance of a 
common leg. The second term describes the coupling due 
to an additional junction with critical current /c,ij inserted 
in the common leg between both qubits. The many-qubit 
system described by Eq. [T] can be characterized by the set of 
eigenstates and eigenenergies /i = 1, .., 2^, which can 
be obtained from the solution of the equation: i^g = E^\ii). 
As an example the first two levels Eq and Ei are shown in 
Fig. [2] for a two qubit system with anti-ferromagnetic coupling 
(Ji2 > 0). 

in. Qubit readout using a low-frequency 

RESONATOR 

For the readout of the qubit system we use a high-quality LC 
resonator (tank-circuit) consisting of a superconducting coil 
with inductance Lt in parallel with a capacitance Ct [14]. 
In order to analyze the response of such a tank to the qubits 
we have to consider the complete Hamiltonian of the qubits 
system and readout resonator: H = Hq -\- Ht + i^int , where 
Ht is the Hamiltonian of the driven tank circuit and i^int the 
interaction term between tank and qubits. In principle, one 
should also add a term containing the dissipative interaction 
with the environment, which would result in line broadening, 
F^j,, of the many-qubit spectrum. However, the measurements 
are performed at a resonator frequency cot which is much 
less than the energy spacings, co^j^ = E^y : ujt <^ \oj^v\^ of 
the coupled qubits. Under such non-resonant conditions the 



effects of dissipation on the measurements are proportional to 
the small ratio V / \uo and can be neglected. 

The driven tank-circuit Hamiltonian can be written in terms 
of the non-commuting current. It and voltage, Vt operators, 
(here [Vt,^t]- = —ifioor^^ with uot = \JLtCt the 
resonance frequency) as: 

CtVt EtIt 

Ht = \ ^T^T ^drive, (3) 

where /drive (^) = ^ac cos(cjt) is the driving current of the tank. 
The interaction term is given by Hi^t — —It^ with ^ the 
magnetic flux through the tank loop created by all qubits, $ = 
Aicri*\ Xi = MiTlp^i is the coupling parameter between 
the i-th qubit and the tank, which is proportional to the mutual 
inductance between qubit and the coil of the tank-cirquit MiT 
and to the magnitude of the persistent current of the qubit. By 
using the Heisenberg equations ihVT = [Vt, H]- and ihiT = 
[It jH]- we obtain the following equation for the averaged 
current and voltage in the LC-circuit: (It) = {Vt)/Lt-, 

(4) 

where the damping of the tank jt — ^t/Qt was added on 
an ad-hoc basis. The flux <l>tot incorporates not only the flux ^ 
which the qubits would apply without their interaction with the 
tank, but also includes the response due to the tank according 
to: 

^tot = m^{^/^) f dt^[m.Hh)]-lT{h)' (5) 

J — oo 

Here we have used the following notations for the 
(anti)commutator of two operators A and B: [A, 5]± = AB^ 
BA^ and the interaction picture. 

Eq.Hlhas the solution {VTit)) = Vt cos(a;t + 6) with 6 the 
phase shift of the voltage relative to the driving current. For the 
case of resonant driving, uj = ujt, the voltage-current phase 
shift is determined by the sum, over all qubit combinations, 
of the real part of the qubit susceptibility Xij • 

tan6 = -|^^A,A,x-,(^T). (6) 

The many-qubit susceptibility at the frequency of the tank 
cotj^t {E^ — Ejy)/h is calculated in the appendix and 
results in: 

^^. = E|^(^i^^'V)(H^i^'V)- (7) 

Here , denote the equilibrium occupation of the qubit 
energy eigenstates: = e"^^/^^^/(X]^ e"^"/^^^). A 
numerical calculation of the predicted response for a 2 qubit 
system is shown in Fig. [3] By comparing this with the energy 
diagram in Fig. [2] it is seen that the predicted signal is largest 
close to the anticrossings where the ground state Eq has the 
largest curvature due to the superposition of flux states. 
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Fig. 2. The first two energy levels for an AFM-coupled two qubit system. 
Calculated for = 100 nA, /p,2 = 150 nA, Ai = A2 = 100 mK and 
J=300 mK. 
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Fig. 3. — tan O for an AFM-coupled two qubit system. The tank response 
is calculated for the same qubit parameters as used in Fig. |2] The read-out 
parameters are Lt = 80 nH, Mi^t = ^2,T = 50 pH and Qt = 500 with 
an effective temperature T = 70 mK. 



IV. Experimental Implementation 

For our experiments we fabricate the qubit-system of in- 
terest inside a prefabricated superconducting planar coil made 
out of niobium on a Si chip. This coil forms the tank-circuit 
together with an external capacitor [15]. On the coil chip 
there are also Nb lines which allow the application of an 
asymmetric flux bias to the qubits. By applying a dc bias 
current /^t through the tank and currents /51 and (in the 
3 qubit example shown in Fig. [T] |7^) through the lines we can 
separately change the magnetic flux in the qubits. For more 
qubits one should add additional lines. 

As shown above, measuring the shift of the resonance 
frequency allows one to read-out the total susceptibility of 
the qubit system. Therefore, the tank-circuit is driven with a 
current /drive (^) = IacCos{u;t) at a frequency co close to the 
resonance frequency cot of the tank circuit. The variation of 
the resonance frequency, due to the change in the susceptibility 
of the qubit system according to ©, can be measured by 
detecting the phase shift between the ac-current /drive 
applied through the coil and ac-voltage over the tank circuit. 
This is realized by using the setup shown in Fig.[T] The voltage 
across the tank circuit is amplified by a cold HEMT-based 
amplifier followed by a room temperature amplifier. The phase 
shift B is measured using a lock-in amplifier with part of the 
signal which is used to bias the qubit used as a reference. 



In the simplest case - one qubit at low temperature ksT <^ 
{E^ — Eiy)/h, ^tot is proportional (through the mutual in- 
ductance) to the persistent current Iq flowing in the qubit 
system. Since a qubit is in the ground state £^0 the current 
Iq is proportional to . By taking into account Eq. |4] it can 
be easily shown that a shift of the tank resonant frequency 
due the tank-qubit interaction (and, therefore, tan 6 or Xij ) 
is proportional to the ground state curvature Because 
the ground state curvature is maximal near the anticrossing, 
the output signal as a function of the applied magnetic flux 
undergoes a narrow dip [16] at <l>o/2. From the shape of the 
dip the one-qubit Hamiltonian can be reconstructed by making 
use of the following equations [8]: 



tan Omax^oFHW 



/. 



^2V22/3_i' 



(8) 



(9) 



At higher temperatures one can fit the dip directly to the theory 
as given by Eqs. (|6]) and (|7]). 

By properly biasing a multi qubit system one can use this 
single qubit reconstruction procedure for a multi qubit system. 
The bias to such an N qubit system is applied as follows: for 
N — 1 qubits and Jij for any j. Then the 

in Eq. ([T]) for — 1 qubits and also the interaction between 
all qubits can be neglected. Doing so reduces Eq. ([T]) to the 
Hamiltonian of a single qubit. The A^ — 1 qubits far from their 
degeneracy points can be considered as conventional classical 
magnetic moments and their influence on the Nth qubit pro- 
vides an additional bias eAr,eff only. The resulting Hamiltonian 
is Hq = — (eAr,eff + eAr)cri^^ + Aatct^^^ + Const. Therefore, 
from measurements described above, and by making use of 
Eqs. dSJ), ©, the persistent current Ip^N as well as A at can 
be determined. By repeating this procedure for the other qubits 
we can find their parameters. 

In order to complete the determination of the parameters 
of Hq the coupling energies Jij must be obtained. They 
can be obtained from the shift of the qubit dips relative 
to the common degeneracy point. This can be seen in the 
numerically calculated energy levels and qubit response in 
Figs. [2]and[3l For low T and small A^, the locations of the 
dips in tan 0(/i, /2, • • •)' anti-crossings, simply follow 

the separation lines between the different classical states in 
the classical stability diagram. Due to this the state of the 
two-qubit system can be easily reconstructed. Moreover, the 
coupling energy follows directly from this stability diagram. 
For instance, the transition |||) < — > |||) occurs at ei = J. 
Therefore, one can easily show that the coupling energy 

follows from J = (5/gl^o/p,l = ^fq2^oIp,2' 

In Fig. |4] we show some quite recent results from a 2 qubit 
system measured with an improved setup. The qubits in this 
sample are coupled by a common Josephson junction as in 
Ref. [12]. This allows one to design the coupling strength 
over a wider range and than with inductive coupling only. 
By using the procedure described above all parameters of the 
Hamiltonian, Eq. [T] were found from the experimental data. 
The mutual inductances M^t between the coil and qubits were 
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Fig. 4. (a) — tan B (7^2, /bT) for AFM-coupled qubits measured at the 
mixing chamber temperature of 10 mK. (b) Theoretical fit for an effective 
temperature Teff = 30 mK and for the sample parameters Ai = 45, A2 = 
55, J = 550 (all in mK) and Ipi = 92 nA, Ip2 = 84 nA. 




Fig. 5. Energy of the system for various vectors. J12 = J23 = J13 = 0.3 K. 
ei = 0.315 K, 82 = 0.252 K, and £3 = 0.525 K. iTibts) is a global 
minimum, while |iit2T3} is a local one. 



determined from the periodicity of the qubit signal from a scan 
over multiple ^o- The mutual inductances between the bias 
lines and qubits followed from the slope of the qubit lines in 
Fig. [4^. These mutual inductances can be used to calculate 
the energy biases and, therefore, for the prediction of the 
response of the tank from the theory of Sec. Hill This prediction 
is shown shown in Fig. |4]3. The good agreement between 
experimental data and the theoretical prediction confirms the 
systems effective temperature of about 30 mK which also 
followed from the fit. 

V. Adiabatic Quantum Computation 

In this section we review the AQC-algorithm for the solution 
of the MAXCUT problem. In particular we discuss a possible 
demonstration of AQC by making use of a system composed 
of three coupled flux-qubits using a resonant tank circuit for 
the readout. 

Since the Hamiltonian of Eq. ([T]) is similar to an Ising 
Hamiltonian at A. = it encodes the MAXCUT of a 3 vertex 
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Fig. 6. — tan B(/q2, fqs) at fqi = 0.006. The tank's response is calculated 
for Ai = A2 = A3 = 96 mK, J12 = J13 = J23 = 300 mK, Ipi = 
Ip^ = 350 nA, Ip2 = 420 nA, and T = 10 mK. The white dashed lines 
denote the cross-overs between the different classical states. At the white 
dot (0.006, 0.004, 0.01) the MAXCUT problem with the solution |TU2T3> 
is encoded, see Fig. \5\ The solid arrows show the directions in which the 
read-out should be carried out in order to reconstruct this state. 



graph if Ai/Ji <C 1 [8]. The simplest non-trivial case is a 
three-qubit system as shown in Fig. [TJ for example with the 
following qubit parameters: Ai = A2 = A3 = 96 mK, J12 = 
Ji3 = J23 = 300 mK, Ipi = Ip3 = 350 nA and Ip2 = 420 nA. 
If this system is allowed to adiabatically evolve to the qubit 
energy biases si (0.006) = 0.315 K, £2(0.004) = 0.252 K, 
and £3 = 0.525(0.01) K, it will encode a MAXCUT problem 
with a solution given by the state |Tii2T3). This can be seen in 
Fig. [5] where the energies of all flux-states for the Hamiltonian 
of Eq.O with all A^ = 0, Jij and Ei as given above have been 
depicted. This system also exhibits a local minimum |iiT2T3), 
two spin-flips away from the global minimum, this makes it an 
interesting test system for AQC. If a violation of the adiabatic 
evolution criterium occurs, the system can be found in this 
local minimum instead of the global one. This can be due to 
either thermal excitations during the readout or Landau-Zener 
transitions caused by a too high readout speed. Therefore one 
should optimize this readout speed depending on the energy 
gap in order to obtain optimal results. 

The tank circuit response — tan 6(/g2, /gs) predicted by 
Eqs. © and d?]) for a fixed value fqi = 0.006 is shown in 
Fig. [6l The white dashed lines denote the cross-overs between 
the different classical (A^ = 0) states, which are marked in 
the figure by the state vectors |iiT2T3), |TiT2i3), |Tii2i3) 
and 111 1213). The main feature of the tank's response is that 
it restores the classical cross-overs. This fact can be used for 
reading out the state of the qubits at any flux point. In the 
case of Fig. [6l we want to know the qubits' configuration at 
(0.006, 0.004, 0.01), because this point encodes the solution of 
the MAXCUT problem that we are interested in. Therefore, 
if we perform the measurement of the qubits' susceptibility 
in the directions marked by solid arrows in Fig. [6l a peak 
structure will appear either to the left, or the right side of 
(0.006,0.004,0.01) in each flux direction. This allows the 
establishment of the following criterium: if the peak appears 
to the right(left) side of the starting point when we scan the 
flux fqi through the qubit i, the starting point corresponds to 
the classical state li (ti). Therefore, one can easily check 
that (0.006,0.004,0.01) corresponds to the |Tii2T3) state if 
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Fig. 7. (a) micrograph of the sample (b) — tan B(/g2, /gs) at fqi = 0.008, 
measured at Tmix = 10 mK. (c) predicted response forAi=A2 = A3 = 
70 mK, J12 = Ji3 = J23 = 610 mK, Ipi = Ip2 = 115 nA, Ip3 = 125 nA, 
and Teff = 70 mK. (d) Difference between the ground and first excited states. 
For (c) and (d) the black dashed lines denote the cross-overs between the 
different classical states. 



all = (for finding /i one would need another figure e.g. 
the (/qi, fqs) plane at a fixed value of fq2 = 0.004). 

As a first step towards the demonstration of AQC we 
measured the full three-dimensional susceptibility of three 
anti-ferromagnetically coupled qubits [17]. For this purpose 
we fabricated a sample with the layout of Fig. [T] and measured 
it in our setup at a base temperature of 10 mK. Figure [7^ 
shows a micrograph of the sample: three Al persistent current 
qubits are fabricated inside a Nb pancake coil. Two junctions 
in each qubit are nominally 600 x 120 nm^, while a third 
one is ~35% smaller. Each qubit is coupled to the other two 
both magnetically and via shared 120 x 2000 nm^ junctions. 
The flux through the qubits is controlled by direct currents 
through the coil /bx, and two additional lines /bi and /b2 
(the third line was not used during the experiment). The Nb 
coil has an inductance Lt = 134 nH, and together with 
an external capacitance Ct = 470 pF forms a parallel tank 
circuit with ujt/^tt = 20.038 MHz and quality Qt = 700. 
The qubit-coil mutual inductances were extracted from the ^0 
periodicity of the ac-susceptibility of the individual qubits as 
MqiT ^ 45.8 pH and M^st ^ 46.6 pH, M^st ^ 45.8 pH. 
The mutual inductances between all qubits and biasing wires 
were found from two scans {hr^ hi)-, {hr-, ^2), see Ref. [17]. 

Fig. (TJ) shows - tan6(/g2, /gs) at fqi = 0.008, while 
Fig. [7]: is the theoretical prediction for Ai = A2 = A3 = 
70 mK, J12 = Ji3 = J23 = 610 mK, /^i = 1^2 = 115 nA, 
IpS = 125 nA, and T = 70 mK. The experimental and 



theoretical data are found to be in good agreement in the full 
three-dimensional flux space (other data for fqi = 0.005, 0, — 
0.005, —0.008 can be found in [17]). However, for this sample 
not all classical cross-overs could be reconstructed. The narrow 
lines in the upper right quadrant of Fig. |7]3C are absent. There 
are two main reasons for this: firstly the difference between the 
ground and first excited state (35 mK) along these measured 
lines is small (see FigjTjl) in comparison with the effective 
temperature T = 70 mK for this measurements. Secondly all 
qubits' persistent currents (or magnetic moments) are similar, 
therefore there is no significant magnetization change along 
the transitions |TiT2i3) < — > |iiT2T3) and |iiT2T3) < — > 
111 1213), thus there is no susceptibility change and no phase 
shift. In order to observe these lines, one could make a sample 
with different persistent currents, increase the gap between 
the ground and first excited states, or decrease the system's 
effective temperature. For the latter approach we already made 
some progress because for our improved experimental setup, 
see FiglH we have reached an effective temperature ~ 30 mK. 

The next steps for AQC will be a improving the read- 
out speed, the demonstration of the real adiabatic evolu- 
tion/computation and demonstration of its efficiency in com- 
parison with the classical simulated annealing. 

VI. Discussion 

This review is concluded by answering some general ques- 
tions that we have been often asked. The first question is 
whether it is possible to realize AQC with coupled classical 
magnetic moments? Or in other words: do we need quantum 
mechanics for an AQC implementation? In the quantum sys- 
tem there is an energy gap between the ground and the first 
exited states, thus the qubit system can evolve always staying 
in the ground state. A system without level- anticrossings 
exhibits many metastable states, so that evolution in the ground 
state is not possible. 

As the qubits stability diagrams, shown in Figs. [3l [6] are the 
same as would be is expected for coupled classical magnetic 
moments, where is the proof that the system described above 
behaves quantum mechanically? In the classical limit the 
persistent current "qubits" have a magnetic hysteresis for 
any finite qubit's inductance. This hysteresis was observed 
by making use of our method [18]. In this work instead of 
hysteresis we observed a narrow dip, which reflects the level 
anticrossing due to quantum tunneling. Therefore, there is a 
nonzero gap between the ground and the first exited state 
which allows AQC. 

The stability of AQC against external noise, to our knowl- 
edge, is not known. On one hand, since we do not use any 
coherent oscillations, AQC should not be very sensitive to the 
environment. On the other hand, external noise can change the 
Hamiltonian itself and, therefore, its ground state. This issue 
requires further study. 

In the present work we discussed our approach for demon- 
strating an adiabatic quantum algorithm by making use of a 
coupled flux-qubit system. By making use of a low-frequency 
resonator inductively coupled to the qubits we can completely 
reconstruct the parameters of a multi-qubit system. This multi- 
qubit ground state anticrossings read-out can be currently 
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performed at an effective temperature of 30 mK. The ex- 
perimental data are found to be in complete agreement with 
quantum mechanical predictions. We have reconstructed the 
susceptibility of three coupled flux qubits in full parameter 
space. The next steps towards AQC will be improvement of 
the read-out speed and the demonstration of adiabatic quantum 
computation and its efficiency. 

Appendix 
Susceptibility calculation 

Starting from the first order approximation of the total flux 
given in the main text (eq. [5]) we can write out the commutator 
as a sum of the qubit response functions (^ij : 

((V^)[$(t),^(ti)]_)^(t-ti) = ^A,A,(^,,(t,ti). (10) 
with 

^ij{tM) = {{i/n)[a^\t),a'^^\h)].)e{t-U), (11) 

where we introduced the Heaviside step function: 0(t — ti) = 
1, if t > ti, and 0{t - h) = 0, when t < h. 

For small values of the tank current and voltage the deriva- 
tive of the total qubit flux is given by the expression 

j^i^totit)) = / dtiMt,ti){VT){h). (12) 

Putting this into Eq. results in a shift of the resonant fre- 
quency of the tank which is proportional to the susceptibility 
and is just the Fourier transform of the response function: 



(13) 



This shift in the resonance frequency causes the phase shift 
when the tank is driven with its unloaded resonance frequency 
ujt according to Eq. ©. 

In order to calculate the linear response functions (pij{t, ti) 
(Eq. (fTTt ) we have to find the product of two projections of 
Pauli matrices taken at times t and ti 

{ai^\t)ai^\h)) = J2pM'^i'\t)<^i'HhM, (14) 

averaged over the quantum-mechanical states and the 

equilibrium distribution = e~^^^^^^ /{^^ ^-E^/kBTy 
As a result we obtain the expression 



(Pij{t,ti) = (pij{t - ti) = (i/n)o(t - ti) 



(0 |//\p*^/3m(^-^i) 



(15) 



for the linear response functions (pij{t^ti) STSi . which give 
after Fourier transformation: 



hu^E^-Ey^ zO 



(16) 



By taking into account that the resonant frequency of the tank 
is much smaller than the level spacing and, therefore, can be 
neglected Eq. d?]) of the main text is generated. 
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